Extensive air showers with TeV-scale quantum gravity 



Luis Anchordoqui, Haim Goldberg, Thomas McCauley, Thomas Paul, Stephen Reucroft, and John Swain 
Department of Physics, Northeastern University, Boston, MA 02115, USA 

One of the possible consequences of the existence of extra degrees of freedom beyond the electroweak 
scale is the increase of neutrino- nucleon cross sections (o~ u n) beyond Standard Model predictions. 
At ultra-high energies this may allow the existence of neutrino-initiated extensive air showers. In 
this paper, we examine the most relevant observables of such showers. Our analysis indicates that 
the future Pierre Auger Observatory could be potentially powerful in probing models with large 
compact dimensions. 
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Recently, it has become evident that a promising route 
to reconcile high energy particle physics and gravity is 
to modify the nature of gravitational interactions at dis- 
tances shorter than a millimeter. Such a modification can 
be most simply achieved by introducing extra dimensions 
in the sub- millimeter range |jj . In this approach the fun- 
damental scale of gravity M* can be lowered all the way 
to O (TeV), and the observed Planck scale turns out 
to be just an effective scale valid for energies below the 
mass of Kaluza-Klein (KK) excitations. Clearly, while 
the gravitational force has not been directly measured 
below the millimeter range, Standard Model (SM) inter- 
actions have been fairly well investigated at this scale; 
so if large extra dimensions really exist, one needs some 
mechanism to prevent SM particles from feeling those 
extra dimensions. Remarkably, there are several possi- 
bilities to confine SM fields (and even gravity) to a 4 
dimensional subspace (referred to as a 3-brane) within 
the (4 + n) dimensional spacetime || . The provocative 
new features of this scenario have sparked a flurry of ac- 
tivity to assess its experimental validity. A brief resume 
of current theoretical work devoted to higher dimensional 
models includes topics addressing fundamental issues of 
phenomenology p), cosmology Q, astrophysics ||, and 
gravity Moreover, an intense effort to find signatures 
of extra-dimensions in collider data is currently underway 

0-. 

Since 1966, a handful of extensive air showers have 
been observed corresponding to what seem to be sin- 
gle particles carrying over 10 20 eV Q. This, in itself, 
is remarkable, as it is difficult or even impossible to ex- 
plain how such energies can be attained by conventional 
acceleration mechanisms Deepening the mystery, it 
was pointed out by Greisen, Zatsepin and Kuz'min |lC| ] 
(GZK) that extremely high energy (> 10 20 eV) cosmic 
rays, if nucleons and/or nuclei, would lose energy rapidly 
through interactions with the cosmic microwave back- 
ground (CMB). This leads to the so-called GZK cutoff, 
which limits the propagation distance of these particles to 
roughly 50 Mpc. The difficulty in constructing nearby as- 
trophysical sources that could accelerate particles to such 
high energies led to the belief that beyond roughly 10 20 
eV, no cosmic rays would be detected. Adding to the puz- 
zle, the arrival directions of these events are distributed 
widely over the sky, with no plausible optical counter- 



parts (such as sources in the Galactic plane or in the Lo- 
cal Supercluster). Furthermore, the "super-GZK" data 
are consistent with an isotropic distribution of sources 
in sharp contrast to the anisotropic distribution of light 
within 50 Mpc from Earth Jll]]. In conclusion, the cur- 
rent picture is very unclear. Thus, it is reasonable to 
consider whether new physics could be at play. 

Of particular interest here, the extraordinarily high 
center-of-mass (cm.) energies achieved at the top of the 
atmosphere are well above those necessary to excite the 
hypothetical KK modes which would reflect a change in 
spacetime dimensionality |l^ ]. Hence, a detailed analysis 
of extensive cosmic ray showers, taking into account this 
departure from previous fundamental particle theory, is 
worthwhile [ p^[ . 

Interestingly enough, if gravity becomes strong at ener- 
gies of a few TeV, virtual graviton exchange can produce 
relatively large effects on the high energy scattering cross 
section, drastically changing the neutrino-nucleon inter- 
action fl4|| . Neutrinos can propagate through the CMB 
essentially uninhibited, breaking the GZK barrier Jl5| . 
Unfortunately, within the SM scenario a neutrino inci- 
dent vertically on the atmosphere would pass through 
it unihibited as well, never initiating an extensive air 
shower. It was already noted that within the extra di- 
mensional framework, the neutrino nucleon cross section 
can approach typical hadronic values at cm. energies 
s > 400 TeV, allowing earlier development of a verti- 
cal neutrino induced shower p6HT8||, One may wonder 
whether the growth of the cross section carries with it 
observable deviations from SM predictions. Consistency 
with current experimental data requires |19| , 
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and this bound certainly does not challenge the neutrinos 
acquiring a hadronic-scale cross section. 

A complete theory of massive KK graviton modes is 
not yet available, making it impossible to know the ex- 
act cross section at asymptotic energies. Any air shower 
analysis would thus depend on reliable guesswork, sup- 
plemented with generally acceptable theoretical princi- 
ples such as duality, unitarity, Regge behavior and par- 
ton structure. A simple Born approximation to the elas- 
tic z^-parton cross section |17 (which underlies the total 
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i/-proton cross section) leads, without modification, to 
a vp ~ s 2 ■ Unmodified, this behavior by itself eventually 
violates unitarity. This may be seen either by examining 
the partial waves of this amplitude, or by noting the high 
energy Regge behavior of an amplitude with exchange of 
the graviton spin-2 Regge pole: with intercept a(0) = 2, 
the elastic cross section 



dcrei \A R (s,t)\ 2 
dt ~ s 2 
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whereas 



Im[A fl (0)] o(0) . 
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so that eventually a e \ > a tot- Eikonal unitarization 
schemes modify these behaviors: in the case of the tree 
amplitudes |14| the resulting (unitarized) cross section 
a l°p ~ s ' whereas for the single Regge pole exchange 
amplitude, cr£,°' ~ In 2 (s/ s o) p0[ |. However, the Regge 
picture of graviton exchange is not yet entirely estab- 
lished: both the (apparently) increasing dominance as- 
sumed by successive Regge cuts due to multiple Regge 
pole exchange 1 14 j2^| , as well as the presence of the zero 
mass graviton can introduce considerable uncertainty in 
the eventual energy behavior of the cross section. Here- 
after, we work within the unitarization framework of Ref. 
and adopt as our cross section |22[ 
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To simulate the consequences of this for ^-induced air 
showers, we assume that the increase in the cross section 
is driven by the production of minijets j2^], and we adopt 
the SIBYLL package to model the fragmentation region at 
ultra high energies p4j |. In other words, the probability 
distribution for obtaining TV jet pairs (with P^ ot > P™ m , 
where P™ m is a sharp threshold on the transverse mo- 
mentum above which soft interactions are neglected) in 
a collision at energy ^/s is computed regarding ^-nucleon 
scattering as a diffractive shadow scattering associated 
with inelastic processes p5[ . Particle production comes 
after the fragmentation of hypothetical colorless parton- 
parton chains mimicking that of SIBYLL hadron-hadron 
scattering. The reader should keep in mind the crude- 
ness of this approximation. However, the imposed cut- 
off on the soft processes ensures that the inelasticity in 
any ncutrino-nuclcon collision is not much larger than 
y ~ 0.15 pOj ], justifying the use of the SIBYLL package. 
As we discuss below, most of the expected qualitative 
features in the shower can be quite well reproduced. The 
algorithms of aires (version 2.1.1) |^6| are slightly mod- 
ified so as to track the particles in the atmosphere. In 
particular, Eq. (Q) is translated into the neutrino mean 
free path 
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FIG. 1. Longitudinal development of neutrino and proton 
showers for different primary energies and primary zenith an- 
gle 43. 9o. The error bars indicate the standard fluctuations 
of the means. 



via the standard 8 parameter function used in aires, 



1 + Pi u + P 3 u 2 + Pj u 3 

V ~ 1 1 4- P. 7/ . + Pa 7/2 -L P, 7/3 4- f 



gem- 2 . (6) 



K = 



(5) 



1 + P 5 u + P 6 u 2 + P T u 3 + P 8 W 

Here TO a ir [g] is the mass of an average atom of air, and 
u = In E [GeV]. The coefficients Pi arc listed in Table I 
for different values M». 

Several sets of neutrinos were injected at 100 km above 
sea level. The sample was distributed in the energy range 
of 10 20 eV up to 10 21 eV, and was uniformly spread in 
the interval of 0° to 60° zenith angle at the top of the 
atmosphere. All shower particles with energies above the 
following thresholds were tracked: 750 keV for gammas, 
900 keV for electrons and positrons, 10 MeV for muons, 
60 MeV for mesons and 120 MeV for nucleons. The re- 
sults of these simulations were processed with the help of 
the aires analysis package. 

Figure 1 shows the total number of charged particles 
versus atmospheric depth averaged over 25 showers for 
the case of a 300 EeV neutrino at M*=l TeV. For com- 
parison, proton-induced showers at 60 and 90 EeV are 
shown on the same figure. As showers initiated by neutri- 
nos typically start later than proton-induced showers, the 
longitudinal development tends to level off after reaching 
a maximum, in contrast to a standard air shower which 
decreases more rapidly after reaching a maximum. The 
number of charged particles produced in the cascade de- 
pends on the amount of energy deposited in the atmo- 
sphere by the primary. Neutrinos at the energy and mass 
scale shown in the figure typically suffer 2 interactions in 
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FIG. 2. Lateral distributions of vertical 300 EeV neu- 
trino-induced showers (triangles), 60 EeV proton-induced 
showers (circles) , and 60 EeV iron-induced showers (squares) . 
The error bars indicate the RMS fluctuations. 



the atmosphere; any energy remaining after this is unde- 
tected. By comparing the neutrino-induced showers to 
the proton-induced showers shown in the figure, one can 
roughly estimate the inelasticity to be 0.1 < y < 0.15. 
This is consistent with the estimates of reference pof J*| 

Figure 2 shows the lateral distributions for vertical 
showers produced by 300 EeV neutrinos, 60 EeV pro- 
tons, and iron nuclei of 60 EeV. At 50 m from the core, 
the ratio of the number of charged particles in the neu- 
trino shower to that in the proton shower is ~ 2, whereas 
it is sa 1.5 in ^/ 56 Fe showers. At about 1 km from the 
core these ratios reduce to w 1.1 and « 0.7, respectively. 
This is significant since experiments which rely on sur- 
face detectors to determine shower parameters typically 
use samples taken on the order of 1 km from the core, 
and thus would not be able to easily distinguish between 
these particle species. 

Figure 3 shows the radial dependence of the mean ar- 
rival time of muons for showers initiated by 300 EeV neu- 
trinos and 60 EeV protons. It can be readily seen from 



*It is important to stress that the maximum number of 
charged particles produced in a proton-induced shower does 
not depend on the hadronic interaction model J27|, making 
the present estimate on the inelasticity quite reliable. 
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FIG. 3. Arrival times for charged particles in vertical 
300 EeV neutrino and 60 EeV proton showers normalized at 
50 m from the shower core. The error bars indicate the RMS 
fluctuations. 



the comparison that the proton-induced showers exhibit 
larger fluctuations than the neutrino-induced showers. 
Besides, each profile presents a well defined slope that 
characterizes the shower front and comprises a signature 
of the primary species. In particular, a neutrino interacts 
in the atmosphere only once or twice, and consequently 
the muons reach the ground with a relatively short time 
delay. 

The simulated neutrino showers discussed so far de- 
posit far less energy in the atmosphere than the most 
energetic of the observed cosmic ray events. A natu- 
ral question is then what the shower profile would look 
like for a neutrino whose energy and mean free path are 
such that it would deposit roughly the same energy as 
observed in the highest energy event . 

At this stage, it is important to point out that within 
the SM framework neutrinos are produced at extremely 
high energies, typically by the weak decay of pions or 
other hadrons. Thus, one needs protons to be accelerated 
to energies a few orders of magnitude even higher. In 
scenarii involving precocious unification [ p9| , there may 
be alternatives to decay chains for producing super-GZK 
neutrinos at the source. 

Figure 4 shows the longitudinal development of a 
900 EeV neutrino-induced shower with a fundamental 
mass scale M* = 1.3 TeV. We stress that such a scale is 
above the lower bound for M* derived from the expected 
flux of neutrinos and current non-observation of horizon- 
tal air showers [^2|. The total energy deposited in the 
atmosphere (after 2 interactions) is of the same order as 
the Fly's Eye event, but the shower maximum occurs, as 
expected, significantly later. 

In summary, it has been proposed fl6|-|l8|| that the 
GZK cutoff can be skirted if the progenitors of the most 
energetic air showers are neutrinos. Under this hypoth- 
esis, the neutrino- nucleon cross section is increased by 
the presence of extra dimensions, allowing the neutri- 
nos to interact in the atmosphere. Simulations indicate 
that neutrino-induced showers at energies of a few hun- 
dred EeV would exhibit signatures distinct from those of 
proton (or nucleus) induced showers that deposit a sim- 
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FIG. 4. The longitudinal development of a 900 EeV neu- 
trino-induced shower is shown together with the experimental 
data reported by Fly's Eye. The error bars in the simulated 
points indicate the standard fluctuations of the means. 



ilar amount of energy in the atmosphere. Similarly, if 
there are neutrinos energetic enough to deposit as much 
energy in the atmosphere as is observed in the highest 
energy events, it appears they too may have unique sig- 
natures. In fact, any physics beyond the standard model 
that increases the neutrino-nucleon cross section should 
affect shower observables like longitudinal profile (mea- 
sured with fluorescence detectors) and ground particle 
distributions (measured with surface detectors). This 
article contains some qualitative discussion of relevant 
observables of neutrino-induced showers. As far as we 
are aware, no showers have been observed which are con- 
sistent with these features. If candidates are eventually 
discovered, of course it will be necessary to carry out a 
much more detailed simulation than the one presented 
here. We note that future hybrid detectors such as the 
Pierre Auger Observatory |5(| will be in an exceptional 
position to search for such phenomena. 

Note added: After this paper was written, it was 
stressed that extremely high energy (300 EeV) neutrinos 
with larger cross section (s 2 rise) can create showers that 
would look like the highest energy event pi] . If this is 
the case, it should also be stressed that neutrinos of a few 
tens of EeV could induce vertical air showers with very 
distinctive profiles. In Fig. 5 we show the longitudinal 
development of showers initiated by neutrinos of E = 5 x 



10 19 eV.[] For comparison we also show showers induced 
by gamma rays and protons of E = 5 x 10 18 eV. It is 
easily seen that within this framework a 50 EeV neutrino 
shower presents its own signature J32|. 

The question of whether the interaction cross section 
of neutrinos with matter could be greatly enhanced (via 
massive spin-2 exchange) at high energies, is yet unde- 
cided. Observation of deeply penetrating showers with 
5 x 10 18 eV deposited in the atmosphere would give an 
experimental and definite answer to this question. As 
an immediate spinoff, we have the converse fact, i.e., 
that if there were no possible candidate which could 
be associated with a neutrino shower, then it should 
be understood as a serious objection to the hypothesis 
of neutrinos as progenitors of the "super-GZK" events. 
We strongly recommend that the Fly's Eye data be re- 
analyzed searching for evidence of neutrino showers. 
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TABLE I. Coefficients for mean free path parametrization 
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